Spin-wave theory at constant order parameter 
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We derive the low-temperature properties of spin- 5 quantum Heisenberg magnets from the Gibbs 
free energy G{M) for fixed order parameter Al. Assuming that the low-lying elementary excitations 
of the system are renormalized spin waves, we show that a straightforward 1/S expansion of G{M) 
yields qualitatively correct results for the low-temperature thermodynamics, even in the absence of 
long-range magnetic order. We explicitly calculate the two-loop correction to the susceptibility of 
the ferromagnetic Heisenberg chain and show that it quantitatively modifies the mean-field result. 
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I. INTRODUCTION 

For many years the magnetically ordered state of quan- 
tum Heisenberg magnets has been studied with the help 
of the spin- wave expansion.^ This expansion is usually 
implemented by expressing the components of the spin 
operator Si = {Sf,Sf,Sf) at lattice site i in terms of 
canonical boson operators hi and b], using either the 
Holstein-PrimakofF transformation^ or the Dyson-Maleev 
transformation^i^ For example, the Dyson-Maleev trans- 
formation for a spin-i* ferromagnet is 



S+ = (25') 1 _6T^^/25 



s; 



S^ = S-blb,, 



(la) 

(lb) 
(Ic) 



where Sf^ = Sf ±iSf, and the number of bosons at a lat- 
tice site may not exceed 2S in order to faithfully represent 
the 25* -|- 1 eigenstates of S'f . If long-range magnetic order 
is present in the z direction and the number of bosons is 
small compared to 25*, the two-body bosonic interaction 
arising from the term in Eq. Hla|l involving three bosonic 
operators can be treated within conventional many-body 
perturbation theory, controlled by the formal small pa- 
rameter 1/5'. Due to the maximum-occupancy constraint 
the effective bosonic Hamiltonian implicitly contains pro- 
jection operators that eliminate the unphysical part of 
the bosonic Hilbert space. Usually, these projection op- 
erators are simply ignored. For Heisenberg ferromag- 
nets in three dimensions it was shown by Dyson^ that 
in the thermodynamic limit the low-temperature ther- 
modynamics can indeed be obtained without taking into 
account the so-called kinematical interactions associated 
with these projection operators. 

While in the 1960s and 1970s ordered magnets have 
been intensely studied, in recent years the center of at- 
tention has shifted to low-dimensional magnets without 
broken symmetries. In this case the conventional spin- 
wave approach described above is not applicable, because 
it relies on the existence of long-range magnetic order. 
Nevertheless, in many magnetic materials the elemen- 
tary excitations still resemble the spin waves of an or- 



dered magnet. For example, in two-dimensional quan- 
tum Heisenberg ferromagnets^ and antifcrromagncts^ at 
low but finite temperatures, where the order parameter 
correlation length ^ is exponentially large, spin waves 
with wave vectors |fc| 3> C ^ are well-defined elemen- 
tary excitations.'' Other examples for systems where the 
low-energy physics is dominated by elementary excita- 
tions of the spin-wave type are Haldane-gap antiferro- 
magnets (i.e. one-dimensional Heisenberg antiferromag- 
nets with integer spin 5') and one-dimensional Heisenberg 
ferromagnets with arbitrary spin. 

To study the low-temperature properties of these 
systems, several methods have been proposed. The 
Schwinger-boson mean-field theory of Arovas and 
Auerbachr- is perhaps aesthetically most appealing. How- 
ever, going beyond the mean-field approximation within 
the Schwinger-boson approach has turned out to be quite 
difficult At the mean- field level the modified spin- wave 
theory (MSWT) proposed by TakahashiiS, is an alterna- 
tive to the Schwinger-boson approach. MSWT yields re- 
sults that agree with the predictions of Schwinger-boson 
mean-field theory up to numerical prefactors. This is not 
surprising, because both approaches are in fact equiva- 
lent to a one-loop renormalization group calculation's^ 
Recently, Takahashi's MSWT has also been applied to 
more complex problems, such as frustratedii or disor- 
dered magnets, 12 or magnetic molecular clustersii^ How- 
ever, the MSWT has shortcomings: (i) it is very difficult 
to systematically calculate corrections due to interactions 
between spin waves within MSWT and (ii) the absence 
of long-range magnetic order is not obtained as a result, 
i.e., the magnetization is set to zero by hand; this leads 
to ambiguity in the choice of the constraint if the MSWT 
is applied to systems with more complicated magnetic or- 
der, such as ferrimagnets»ii In this work we shall show 
that these problems can be resolved within the conven- 
tional spin-wave approach simply by performing the cal- 
culation at constant order parameter. 

This paper is organized as follows. In Sec. ^1 we dis- 
cuss the calculation of thermodynamic observables at 
constant order parameter. In Sec. lIIII this approach is ap- 
plied to the Heisenberg ferromagnet in D — 1,2,3 dimen- 
sions within linear spin-wave theory. Hartree-Fock and 
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two-loop corrections are obtained for the one-dimensional 
case in Sec. IIVI The work is summarized in Sec. ]V\ 



II. THERMODYNAMICS AT CONSTANT 
ORDER PARAMETER 

In this section we discuss the calculation of thermo- 
dynamic observables at constant order parameter. Al- 
though this approach is applicable to a variety of corre- 
lated systems with order parameter M and correspond- 
ing conjugate field h, here we will focus on the spin-S* 
Heisenberg ferromagnet with zero-field Hamiltonian 



H 



(2) 



where the sum is over all nearest-neighbor pairs of a D- 
dimensional hypercubic lattice with N sites, and J > is 
the exchange coupling. In this case the order parameter is 
simply the total magnetization, given by M — ^i: 
and h is the homogeneous magnetic field (in suitable 
units). Applications to antiferromagnets or more com- 
plicated magnetic systems are straightforward. 

Let us first recall some elementary thermodynamics. 
For fixed field h and temperature T, thermodynamic ob- 
servables can be obtained from the Helmholtz free energy 
(setting the Boltzmann constant to unity) 



(3) 



where the dependence on T is suppressed for brevity. 
Given F{h), the magnetization is obtained as 



M{h) 



dFjh) 
dh 



(4) 



Alternatively, we may choose to fix the magnetization 
and adjust the magnetic field appropriately. The corre- 
sponding thermodynamic potential is the Gibbs free en- 
ergy G{M), which is related to the Helmholtz free energy 
via a Legendre transformation^ 



G{M) = h{M)M + F{h{M)) 



= -TlnTr 



-[H-h{M)(M-M)]/T 



(5) 



where the function h{M) is obtained from Eq. |@}. From 
G{M) we obtain the equation of state in the form h = 
h{M) via 



h{M) 



dG{M) 
dM 



(6) 



which shows that the equilibrium magnetization for van- 
ishing field is an extremum of G{M). If the system has a 
finite spontaneous magnetization Mq = lim/j^Q+ M{h), 
then the generic expected behavior of G{M) is, for M > 
Mo, 



G{M) = G(Mo) 



(M - Mof 
2X 



+ 0[{M-Mof], (7) 



while for \M\ < Mq the Gibbs free energy has the con- 
stant value G{Mq); see, for example, Ref. 16. Here 



X = 



dh{M) 



dM 



d^G{M) 



dAP 



(8) 



Mo 



is the inverse longitudinal order parameter susceptibility 
for vanishing external field. These expressions are also 
valid in the absence of spontaneous symmetry breaking, 
where Mo = 0. Note that, in general, G(M) = G{-M), 
because the spectrum of M is symmetric with respect to 
the origin. 

The parameter h{M) in Eq. © can be viewed as a 
Lagrange multiplier that enforces the condition of con- 
stant magnetization. The zero-temperature version of 
the method outlined above has been used previously by 
Georges and Yedidiaii to study spontaneous symmetry 
breaking in the ground state of the Hubbard model. Note 
that in the limit T ^ Eq. (jsj can be written asi^ 
G(M) = (0|G(Af)|0), where |0) is the ground state of the 
"free-energy operator" G(M) = H - h{M)[M - M]. As 
shown in Ref. llTl the expansion at constant order param- 
eter is advantageous for the calculation of corrections to 
the mean-field approximation. In the following section 
we show that for low-dimensional Heisenberg magnets 
without long-range order this method yields reasonable 
results even at the level of linear spin-wave theory; the 
leading fluctuation corrections in Z) = 1 are then calcu- 
lated in Sec.HVI 



III. LINEAR SPIN WAVES AT CONSTANT 
ORDER PARAMETER 

We now calculate G(M) within linear spin- wave the- 
ory, i.e., to leading order in l/S*, assuming that the low- 
lying elementary excitations of the system are renormal- 
ized spin waves. In this approximation the square brack- 
ets in Eq. Ijlal) are simply replaced by unity, so that the 
Heisenberg Hamiltonian (0) becomes 



Ho = -DJNS^ + J2 ^kbik 

k 

where = 2DJS{\ — 7fc), with 

D 

7fc = -D"^ ^ cos(fc • a^) . 



(9) 



(10) 



For simplicity, we impose periodic boundary conditions 
on a hypercubic lattice with primitive lattice vectors 
and lattice spacing a = |a;^|. The momentum sum is 
over the first Brillouin zone and bk is the lattice Fourier 
transform of bi. The corresponding free energy is 



Fo {h) = -DJNS'^ - hNS + tJ2^^[ 

k 



1 -e 



-{ck + h)/T 



(11) 
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From Eq. we then obtain the usual spin-wave result 
for the magnetization 



M{h) = NS -Y^[e 



{ih + h)/T _ 



(12) 



A. Three-dimensional ferromagnet 

It is instructive to begin with linear spin-wave theory 
for the three-dimensional Heisenberg model. In the ther- 
modynamic limit we obtain for the magnetization per site 
m = M/N to leading order in t = T/JS and v = h/T 



m{h) = S 



87r3/2 ^47r 



(13) 



where C,{z) is the zeta function. Setting /i = we re- 
cover the well-known Bloch T^/^ law for the sponta- 
neous magnetization per site, mo = \\mh^o+ m{h) , in 
the ordered state of the Heisenberg ferromagnet. Tak- 
ing the derivative of Eq. H13|) with respect to /i, we see 
that the susceptibility x = dM/dh [Eq. ©] diverges for 
/i ^ as This diver gence of the uniform longi- 

tudinal susceptibility of a three-dimensional Heisenberg 
magnet in the ordered state is not widely appreciated, 
although it was noticed a long time ago2i^ and has been 
confirmed by renormalization group calculations for the 
classical Heisenberg ferromagnefei2i2£ and perturbative 
calculations for the corresponding quantum model^iiSS 
Due to this divergence, the Gibbs free energy G{M) of 
the Heisenberg ferromagnet in D = 3 does not have the 
generic form Q. Instead the linear spin- wave result for 
G{M) is, for m> mo, 



GojM) 
NT 



T 



C(|) 

87r3/2 



167r 



+ 0[im - ma)^] . 



(14) 



From Eq. (|13|l we note that h^^^ cx (m — rno), so that 
we cannot solve for h as a function of m unless m > 
niQ. In light of the above general discussion this is not 
surprising, because for \m\ < mo the Gibbs free energy is 
constantii^ The behavior of Go{M) as a function of M is 
shown in Fig. ^ The leading m dependence of Eq. H14() 
is proportional to (m — mo)^, which can be traced to the 
fact that the inverse susceptibility vanishes. By contrast, 
in D > 4 the uniform longitudinal susceptibility of the 
Heisenberg ferromagnet is finite, ^^'^"^ so that in this case 
the Gibbs free energy has indeed the generic form O. 



B. One-dimensional ferromagnet 

Let us consider now the one-dimensional case, where 
we know that the Heisenberg ferromagnet does not have 
any long-range order at any finite temperature T. In this 



Go(M) 



-m{h=0) m(h=0) 
m 

FIG. 1: Gibbs free energy Go(M) of the three-dimensional 
Heisenberg ferromagnet within linear spin-wave theory. Due 
to the divergent longitudinal susceptibility in D = 3, the 
Gibbs free energy grows cubically for \m\ slightly above mo, 
see Eq. HM . 



case the linear spin-wave theory result for the magneti- 
zation per site is 

(15) 

where again t = T/JS and v — h/T. This expression 
predicts a divergent magnetization and susceptibility for 
/i — !■ 0. However, we can obtain a perfectly finite result 
for the susceptibility at constant magnetization [Eq. 0]. 
Solving Eq. (|15ll for has a function of M ~ Nm we obtain 



h{M) 



J.2 



4JS[S - 



Oil 



(16) 



According to Eq. (jSJ this implies for the inverse suscep- 
tibility 



rp2 



X 



2NJS[S-m-^Vt] 



(17) 



Anticipating that in one dimension m = 0, we obtain for 
the susceptibility per site at low temperatures 



N 



2JS^ 



1_ 1.^(5; 



S2^ 



Vi + 0{t) 



(18) 



This expression agrees exactljtS^ with the prediction 
of the MSWT advanced by Takahashi^ who argued 
that Eq. H18|l is indeed the correct asymptotic low- 
temperature behavior of the susceptibility for arbitrary 
S. For 5=1/2 the nearest-neighbor Heisenberg chain is 
exactly solvable via Bethe ansatz?^ so that in this case 
one can obtain an independent check of Eq. (|18|l . In- 
deed, from a numerical analysis of the Bethe-ansatz in- 
tegral equations^'* Takahashi^'^ found perfect agreement 
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Go(M) 



Go(0) 




FIG. 2: Solid line: Gibbs free energy Go(M) of the one- 
dimensional Heisenberg ferromagnet within linear spin-wave 
theory, see Eq. 11911 . The cusp at m = is an unphysical 
artefact of the spin-wave expansion, which is related to the 
neglect of the kinematical interaction between the spin waves. 
The dashed line is the subtracted Gibbs free energy Go{M) — 
Go{M) - h{0)\M\. 



with Eq. H18|l for S = 1/2, which is remarkable because 
a priori linear spin-wave theory is only expected to be 
accurate in the ordered state and for large S. We shall 
further comment on this agreement below. 

Within linear spin-wave theory, the Gibbs free energy 
per site is given by 



Go(M) 
NT 



T 
1 

' iVT 



c(i) 



h{Q)\M 



Vt 



0{\M\ 



where x is given in Eq. (|18|l and 



h{Q) 



4JS'3 



0(r5/2) ^ 



(19) 



(20) 



see Eq. (|16f) . In writing Eq. H19|l we have used the fact 
that our spin- wave calculation yields Gq(M) only for 
M > and that the exact G{M) is an even function 
of M . Note that G'o(M) assumes a minimum at M = 0, 
indicating the absence of long-range order. However, as 
shown in Fig. |21 linear spin- wave theory predicts an un- 
physical cusp in the Gibbs free energy at M = 0. The 
finite slope /i(0) = dGaj dM\M=o+ can be identified with 
the variational parameter — /i introduced by Takahashi,^*^ 
which in his calculation plays the role of a chemical po- 
tential for the Dyson-Maleev bosons, enforcing the con- 
dition of zero magnetization. On the other hand, it is 
physically clear that for T > any finite value of the ex- 
ternal field will always be accompanied with a finite mag- 
netization, so that an exact calculation of G{M) should 
yield limM^o^(-^) = 0. Therefore we expect that the 
exact Gibbs free energy in one dimension has the form 



given in Eq. ITJ. The cusp of the spin- wave result for 
the Gibbs free energy is related to the fact that in our 
simple spin-wave calculation we have ignored the kine- 
matical interaction between the spin waves which arises 
from the maximum-occupancy constraint. Fortunately, 
this cusp is irrelevant for the calculation of the zero-field 
thermodynamics, which can be also obtained from the 
subtracted Gibbs free energy 



Go(M) =Go(M)-/i(0)|M|, 



(21) 



see Fig. |21 Note that Gq{M) has the generic behavior 
given in Eq. |7J, with the susceptibility given by Eq. I|18|) . 



C. Two-dimensional ferromagnet 

For completeness we now discuss the case D = 2, where 
the spontaneous magnetization of the Heisenberg ferro- 
magnet is zero at any finite temperature T. The result 
of linear spin- wave theory for the magnetization is (with 
t = T/JS, V = h/T) 

^h) ^S-l-[-\nv + ^ + ^t + 0(t\ «2)] , (22) 

which again diverges for ft ^ 0. The function /i(m) is 
obtained as 



ft(m) = re4^(^-")/*[l + 0(t)]. 



and the result for the susceptibility at m = is 

g47r5/t 



X 



[1 -I- 0{t)] 



(23) 



(24) 



which diverges for T ^ 0. These expressions are anal- 
ogous to Takahashi's resultsiiS Finally, the Gibbs free 
energy takes the form 



Go(Af) _ C(2)^ 4.g/t/ t 
NT All V47r 



TTi \ 

\M + —)+0{t,nv'), 
(25) 



i.e., with a minimum at m = 0, again with an unphys- 
ical cusp; in this sense the situation is rather similar to 
that in D — \. Note, however, that in D — 2 it is known 
that a two-loop calculation is necessary to obtain the cor- 
rect low-temperature asymptotics of the susceptibility.^ 
Although the exponential factor x ^ exp[47r JS'^/T] is 
correctly reproduced by mean-field theory, the two-loop 
correction changes the power of T in the prefactor of 
Eq. (|24l) : the correct low-temperature behavior of the 
susceptibility of the quantum Heisenberg ferromagnet in 
two dimensions is x oc T'^ exp[ATrJS^ /T]. This result is 
not modified if higher-order terms involving more than 
two loops are included..^ 
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IV. BEYOND LINEAR SPIN-WAVE THEORY 



Because fluctuation effects are usually stronger in lower 
dimensions, we expect tliat in one dimension the correc- 
tions to the mean-field result H18|l are even more impor- 
tant than in Z? = 2. Wc now explicitly calculate the 
two-loop correction. Within the Dyson-Maleev formal- 
ism the dynamical spin-wave interactions are contained 
in the following two-body Hamiltonian: 




^ N ^ 



k'^ ,h'^ ,fc2 

xV{k\,k'2MM)hi,bl,hk,K, (26) 



where Sx{k) denotes momentum conservation modulo a 
reciprocal-lattice vector, and the symmetrized interac- 
tion vertex is 



FIG. 3; Feynman diagram describing the leading fluctuation 
correction to the free energy of the ferromagnetic Heisenberg 
model, see Eq. (Onj- The solid arrows denote the Hartree-Fock 
magnon propagators and the squares are the Dyson-Maleev 
vertices. 



y(fci,fei,fe2,fei) = -i 



(27) 



with 7fc defined in Eq. (|10|l . First let us estimate the 
effect of Hi within the self-consistent Hartree-Fock ap- 
proximation. We write our spin-wave Hamiltonian as 
{Hq + 6 Hi ) + {Hi — SHi) and choose the one-body Hamil- 
tonian SHi such that the thermal expectation value of the 
residual interaction Hi — SHi in the ensemble defined by 
the Hartree-Fock Hamiltonian Hq + 5Hi vanishes. We 
obtain 



OF) J 

6H1 = E ^i{k)bibk - -rp E ^('^' ^) 



N 



k,k' 



(28) 



where = [e(^'=+'')/^ — 1]~^ is the thermal occupation 
of the Hartree-Fock magnon states with momentum fc, 
and the Hartree-Fock self-energy is given by 



Si(fc) 



ADJ 

N 



J2vik,k',k',k) 



(29) 



After some standard manipulations we obtain for the 
Helmholtz free energy within self-consistent Hartree-Fock 
approximation 



1 



-(Ek+h)/T 



Fi (h) = -DJNS^ - hNS + T E In [ 

k 

+ DJNS^il- Zf . (30) 

Here Ek — Zek, and the dimensionless renormalization 
factor Z satisfies the self-consistency condition 



1-]4E(1-7.) 



nk ■ 



(31) 



The quantity ZS corresponds to the second variational 
parameter 5" introduced by TakahashiiiS Note that he 
gives a different sign for the last term in Eq. (|30|l . In 
one dimension 1 — Z — 0{T'^) at low temperatures^ 
so that for the calculation of the first two terms in low- 
temperature expansion of thermodynamic observables it 
is sufficient to set Z = 1. We conclude that at the 
Hartree-Fock level the dynamical interaction between 
spin waves does not contribute to the low-temperature 
asymptotics ixi D = \. At this level of approximation 
our theory is equivalent to MSWT. 

Within our approach it is now straightforward to study 
spin-wave interactions beyond the Hartree-Fock approx- 
imation. Therefore we simply expand the Helmholtz free 
energy F{h) to higher order in the interaction and then 
perform a Legendre transformation to obtain the corre- 
sponding Gibbs free energy. We now calculate the first 
fluctuation correction to F{h). The relevant Feynman 
diagram is shown in Fig. 13 In this approximation the 
Helmholtz free energy is F2{h) = Fi(h) + SF2{h), where 
Fi{h) is given in Eq. ll5n|) and 



SF^ih) 



2(—] y W{k,k',q) 



X [{1 + nk)il + nk')nk+qnk'-q 

- UkUk'il + nfc+q)(I -I- Uk'-q)] , (32) 



where 



W{k, k', q) = V{k, k' , k + q,k' q) 

xV{k + q,k' -q,k,k'). (33) 

At low temperatures, we may replace the Dyson-Maleev 
vertex by its long- wavelength limit, which in D dimen- 
sions is given by V{k[, fcj: ^2, ^i) ~ ~{ki-k2)a^ /2D. For 
the rest of this work we shall explicitly set D = 1. Then 
the leading behavior of SF2{h) for small t = T/JS and 
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small V = h/T can be calculated analytically. We obtain 
ior N oo 



5F2{h) 



1 



^3/2 



TN 16 {2SYv^/ 
The resulting equation of state is 

1 A 



— + 0{v"\t^l\). (34) 



U-Jl 

25 25* V V 

3 

1 It ■ 



1 

16 



2SS V 



+ 0{t,t^/\-^/^). (35) 



Comparing this result with the corresponding expres- 
sion obtained within linear spin-wave theory given in 
Eq. H15|l . we see that the two-loop correction gives rise 
to an additional term proportional to the third power 
of {2S)~\t/vy/^. But linear spin-wave theory predicts 
that this parameter is actually close to unity, as is easily 
seen by setting m = in Eq. (|15|l . Hence, the lead- 
ing fluctuation correction to the Hartree-Fock theory is 
not controlled by a small parameter. Note that the extra 
power of S^^ that appears in the two-body part of the 
effective boson Hamiltonian is canceled by the singular h 
dependence of the two-loop correction. If we nevertheless 
truncate the expansion at the two-loop order, we obtain 
from Eq. (|35|l for the leading low-temperature behavior 
of the susceptibility. 



(36) 



with « 1.96, which is slightly smaller than the linear 
spin- wave prediction = 2, and significantly smaller 
than the result C-^ = 3 obtained within Schwinger-boson 
mean-field theory^ We suspect that corrections involving 
more loops will involve higher powers of the parameter 
(2S')~^(i/w)-^/^ in Eq. which give rise to additional 
finite renormalizations of C^. Hence, a numerically accu- 
rate expression for the low-temperature susceptibility of 
a one-dimensional Heisenberg ferromagnet cannot be ob- 
tained from a truncation of the 1/S spin- wave expansion 
at some finite order. Note that quantum Monte Carlo 
simulations for the S* = 1/2 nearest- neighbor Heisenberg 
chain^^ give = 1.58±0.03, supporting the scenario de- 
scribed above. In light of these results it is puzzling that 
from the numerical analysis of the Bethe-ansatz integral 
equations for S = 1/2 TakahashiiSi^ obtained = 2. 
Possibly this is related to difficulties in extracting the 
true asymptotic low-temperature behavior of the suscep- 
tibility from the Bethe-ansatz integral equations pi2i21 



Although at first sight our spin-wave expansion for 
the susceptibility appears to be controlled by the small 
parameter l/S, this parameter is renormalized by the 
infrared singularity of the two-loop correction to the 
mean- field result. This phenomenon is familiar from 
the weak-coupling calculation of the two-loop correc- 
tion to the ground-state energy of the repulsive Hubbard 
model at constant staggered magnetization in one and 
two dimensions Due to infrared singularities inher- 
ent in the loop integrals, those expansions are effectively 
in powers of the Hubbard interaction U multiplied by a 
function of the order parameter; as a consequence the 
two-loop correction to the ground-state energy has the 
same order of magnitude as the Hartree-Fock term. 



V. CONCLUSIONS 

In summary, we have shown that conventional spin- 
wave expansion at constant order parameter is an alter- 
native to Takahashi's modified spin-wave theory, which 
nowadays is one of the most popular mean-field methods 
to study low-dimensional magnets without long-range or- 
der. Due to the conceptual simplicity of our method, we 
can systematically calculate corrections to the mean-field 
approximation using conventional diagrammatic meth- 
ods. We have explicitly calculated the leading fluctuation 
correction to the mean-field result for the susceptibility of 
the ferromagnetic Heisenberg chain. We have found that 
in one dimension the predictions of MSWT are at most 
qualitatively correct, because fiuctuation corrections are 
not controlled by a small parameter. 

Furthermore, compared to MSWT the present ap- 
proach has the conceptual advantage that the intro- 
duction of Lagrange multipliers by hand is not neces- 
sary, since their role is played by external fields instead. 
The difference between these approaches is best visible 
for systems with more complicated order parameters. 
For example, recent work on molecular magnetai^ and 
ferrimagnetsi^ has shown that the formulation of MSWT 
for such sytems is difficult and that the proper choice of 
constraints is not clear a priori. On the other hand, spin- 
wave theory at constant order parameter naturally yields 
the absence of long-range order in low dimensions as a re- 
sult; applications to antiferromagnets and ferrimagnets, 
with vanishing homogeneous and staggered magnetiza- 
tions, are in progress^ 

This work was supported by the DFG via Forscher- 
gruppe FOR 412, Project No. KO 1442/5-1. 



^ See, for example, F. Keffer, in Handhuch der Physik, edited ^ F. J. Dyson, Phys. Rev. 102, 1217 and 1230 (1956). 

by F. Fliigge (Springer- Verlag, Berlin, 1966), Vol. XVIII/2. * S. V. Maleev, Zh. Eksp. Theor. Fiz. 30, 1010 (1957) [Sov. 
^ T. Holstein and H. Primakofl, Phys. Rev. 58, 1098 (1940). Phys. JETP 64, 654 (1958)]. 



7 



^ P. Kopietz and S. Chakravarty, Phys. Rev. B 40, 4858 
(1989). 

® S. Chakravarty, B. I. Halperin, and D. R. Nelson, Phys. 

Rev. B 39, 2344 (1989). 
^ P. Kopietz, Phys. Rev. B 41, 9228 (1990); P. Kopietz and 

G. Castilla, ibid. 43, 11100 (1991). 
® D. P. Arovas and A. Auerbach, Phys. Rev. B 38, 316 

(1988). 

® A. E. Trumper, L. O. Manuel, C. J. Gazza, and H. A. 

Geccatto, Phys. Rev. Lett. 78, 2216 (1997). 
^° M. Takahashi, Prog. Theor. Phys. Suppl. 87, 233 (1986); 

Phys. Rev. Lett. 58, 168 (1987). 
" E. S. Pisanova, N. B. Ivanov, and N. S. Tonchev, Phys. 

Rev. B 65, 212404 (2002). 

X. Wan, K. Yang, and R. N. Bhatt, Phys. Rev. B 66, 
014429 (2002). 

S. Yamamoto and T. Nakanishi, Phys. Rev. Lett. 89, 
157603 (2002). 

S. Yamamoto, T. Fukui, K. Maisinger, and U. SchoUwock, 
J. Phys.: Condens. Matter 10, 11033 (1998). 
See, for example, J. W. Negele and H. Orland, Quantum 
Many-Particle Systems, (Addison- Wesley, Redwood Gity, 
1988), Ghap. 4. 

^® S. K. Ma, Statistical Mechanics, (World Scientific, Singa- 
pore, 1985), Ghap. 27. 



1^ A. Georges and J. S. Yedidia, Phys. Rev. B 43, 3475 
(1991). 

Of course, for a ferromagnet the order parameter and the 
Hamiltonian commute, so that the zero-temperature limit 
is trivial due to the absence of quantum fluctuations. 
E. Brezin and D. J. Wallace, Phys. Rev. B 7, 1967 (1973). 
2" G. Mazenko, Phys. Rev. B 14, 3933 (1976). 

V. G. Vaks, A. L Larkin, and S. A. Pikin, Zh. Eksp. Teor. 
Fiz. 53, 1089 (1967) [Sov. Phys. JETP 26, 647 (1968)]. 
V. P. Kalashnikov and M. 1. Auslender, Physica A 100, 
443 (1980). 

Note that the normalization of the susceptibilty in 
Ref. [13 is different from ours, so that we should identify 
^Takahashi ^ ^^/N. Taking this into account, our Eq. (CHJ 
agrees exactly with the susceptibility obtained by Taka- 
hashi. 

M. Takahashi and M. Yamada, J. Phys. Soc. Jpn 54, 
2808 (1985); M. Yamada and M. Takahashi, ibid. 55, 2024 
(1986). 

M. Kollar, I. Spremo, and P. Kopietz (unpublished). 
P. Kopietz, Phys. Rev. B 40, 5194 (1989). 
P. Schlottmann, Phys. Rev. Lett. 54, 2131 (1985); 56, 1618 
(1986); Phys. Rev. B 33, 4880 (1986). 
^® P. Kopietz, Phys. Rev. B 48, 13789 (1993). 



